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A SIMPLIFIED EQUILIBRIUM HYDROCARBON-AIR 
COMBUSTION GAS MODEL FOR USE IN AIR-BREATHING 
ENGINE CYCLE COMPUTER PROGRAMS 

By Vincent R. Mascitti 
Langley Research Center 

SUMMARY 

This paper presents a simplified hydrocarbon-air combustion gas model, including 
the effects of dissociation, for convenient use in engine cycle computer programs. The 
generalized model reduces to the hydrogen-air system as well as to the dissociating -air 
system. The exclusion of chemical species containing atomic nitrogen allows a consider- 
able simplification of the composition equations. The thermodynamic properties of stoi- 
chiometric combustion of the kerosene-air and hydrogen-air systems are computed with 
the simplified model and compared with those of more comprehensive gas models. In 
addition, the effect of the neglected chemical species on the performance of an idealized 
subsonic combustion ramjet is presented. The simplified gas model has been used to 
define the limiting conditions for solid carbon and ammonia formation for fuel-rich gas 
mixtures. A computer program listing of the calculation procedure for the simplified 
gas model is presented. 


INTRODUCTION 

The current interest in supersonic and hypersonic vehicles has emphasized the need 
for an accurate yet convenient method of calculation of the thermodynamic properties of 
combustion gas mixtures. At high temperatures considerable iteration difficulties arise 
in calculating the chemical composition and thermodynamic properties of combustion gas 
mixtures, difficulties which stem from the phenomenon of internal energy excitation of 
the species coupled with chemical dissociation. 

In the past, methods such as the method presented in reference 1 have been used to 
calculate the thermodynamic properties of combustion gases for air-breathing engine 
cycle computations. The method of reference 1 is based on separate calculations of ther- 
modynamic properties of air and the products of combustion of a stoichiometric fuel-air 
mixture. It is assumed that the properties of the combustion products for any fuel-air 
ratio that is less than stoichiometric may be obtained by linear interpolation between the 



two extreme cases. Results of this method are exact for no dissociation. Application of 
this method to conditions in which dissociation occurs will only give approximate results. 
The deviation from an exact calculation is small for turbojet or timbofan calculations 
since temperatures are relatively low and dissociation is not extensive. However, the 
high operating temperatures of hypersonic engine cycles may result in extreme dissocia- 
tion and cause large errors in calculated engine performance. 

The method of reference 1 does not allow engine cycle calculations for equivalent 
ratios above unity. However, at high Mach number flight, the fuel required to cool aero- 
dynamic surfaces and engine components may force engine operation into this region. 

In short, the method of reference 1 does not provide a sufficiently general basis 
for engine cycle calculations under all conditions of current and future interest. 

Chemical equilibrium theory which governs dissociation phenomenon has been for- 
mulated for many years (for example, ref. 2). However, its application to complex com- 
bustion gas mixtures has proved so tedious and complicated as to preclude direct use in 
engine cycle computer programs. The combustion gas models of references 3 and 4, 
which serve as a basis of comparison for the simplified gas model presented herein, are 
elaborate treatments of the hydrogen-air and kerosene-air systems. 

The generalized hydrogen-air model of reference 3 assumes 12 chemical species 
requiring eight independent chemical reactions leading to eight equilibrium expressions. 
Since the initial proportions of hydrogen and air define four mass balance equations, the 
mathematical solution (12 equations and 12 unknowns) is demonstrated. The reduction of 
the system of equations to one equation and one unknown, as recommended in reference 5, 
is extremely tedious. With the exclusion of the technique of reference 5, solution of this 
system by a single-level iteration is eliminated. The adoption of a two-level iteration or 
one of the methods summarized in reference 6 to solve for gas composition has a great 
effect on the utility of the computation procedure as a subroutine for an engine cycle com- 
puter program. Consequently, the most efficient means of representing the thermody- 
namic properties of these gas models in an engine cycle program is by an elaborate sys- 
tem of fitted curves which sacrifices both computer storage and accuracy. 

The purpose of the simplified gas model presented herein is to incorporate the sub- 
stantial effects of dissociation consistent with convenient use in computerized engine cycle 
calculations. The proposed gas model is simplified by neglecting the formation of species 
containing atomic nitrogen. This assumption enables the solution for chemical composi- 
tion to be obtained with a single -level iteration. Although the original intent of this study 
was directed toward the hydrogen-air system, it was foimd that the model could be gen- 
eralized to any hydrocarbon -air system, as well as to a dissociating -air system, with 
little additional complexity. The computer program, presented in the appendix, can be 
readily incorporated as a subroutine in an ermine cycle program or used alone to generate 
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Mollier diagrams. Inputs to the program are the ratio of carbon atoms to hydrogen 
atoms in the fuel molecule, equivalence ratio, temperature, and pressure. Outputs from 
the program are the mole fractions of the chemical species assumed, enthalpy, entropy, 
and molecular weight. 


SYMBOLS 

A, B coefficients of iteration functions 

D, DijDo,^ 

^ > initial atomic proportions 

D3,D4 J 

^net thrust, Ibf (N) 

f fuel -air ratio 

fg stoichiometric fuel -air ratio 

g acceleration due to gravitational field of earth, 32.174 ft/sec^ (9.807 m/sec2) 

G Gibbs free energy, Btu/lbm (JAg) 

H enthalpy, Btu/lbm (JAg) 

Igp specific impulse, sec 

J mechanical equivalent of heat, 778.20 ft-lbf/Btu (9.991 N-m/J) 

K equilibrium constant in terms of partial pressures; numerical subscripts 

indicate individual reaction 

M mean molecular weight 

n ratio of hydrogen atoms to carbon atoms in fuel molecule 

p pressure, atm 

q dynamic pressure, Ibf/ft^ (N/m2) 
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R 


viniversal gas constant, 1.98588 Btu/mole-°R (8.314 J/mole-°K) 


r^h ratio of carbon atoms to hydrogen atoms in fuel molecule, 1/n 

rgq equivalence ratio, f/fg 

S entropy, Btu/lbm-°R (j/kg-°K) 

T temperature, °R (°K) 

V velocity, ft/sec (m/sec) 

air flow rate, Ibm/sec (kg/sec) 

X mole fraction 

Subscripts: 

oo free-stream conditions 

J Jth species 

j jet 

lim limiting condition for solid carbon or ammonia formation 

s static conditions 

t stagnation conditions 

Superscript: 

o at reference pressure of 1 atmosphere (1.01325 x 10^ N/m^) 

DESCRIPTION OF GAS MODEL 

The temperatures of interest for air-breathing engine application are assumed to be 
below 7000° R (3889° K) and the pressures are assumed to be between 0.001 and 100 atm. 
For these conditions, references 3 and 4 indicate that the formation of nitrogen species. 
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such as N, NH, NH 3 , and NO, occurs in negligible amounts and, therefore, has a very 
small effect on the thermodynamic properties of combustion gas mixtures. The assump- 
tion to neglect these species greatly reduces the complexity of the calculation procedure 
since molecular nitrogen can be considered inert. The results of this assumption on the 
thermodynamic properties of gas mixtures and on the performance of an idealized ram- 
jet are discussed in a subsequent section. 

Figure 1 is a diagram illustrating the chemical species and required chemical 
reactions considered in this gas model. The figure is divided into regions of tempera- 
ture and equivalence ratio. Above T^.^^_Qff, the gas is considered to be dissociating, 
("^cut-off arbitrarily defined limit below which dissociation can be considered 

negligible.) The chemical species assumed for pure air ^r^q = O) are identical to those 
of reference 3 and include atomic nitrogen species. For the combustion gas model 
(req > O), the dissociated nitrogen species ^N, NH, NH3, NO) are not included. It is appar- 
ent that by excluding the species and reaction containing carbon, the model reduces to 
the hydrogen-air system. 

Below there is no dissociation and combustion is complete. Therefore, 

the initial proportions of elements are sufficient to define the gas composition. However, 
for r__ > 1, the gas model of reference 4 indicates that CO may form even below 
Tcut-off- ^ result, the initial proportions of elements are not sufficient to define the 
gas composition and one chemical reaction is required. 

In this region of fuel -rich operation, solid carbon and ammonia may form. The 
limiting pressure for the formation of solid carbon depends upon the relative proportions 
of carbon dioxide CO 2 and carbon monoxide CO, equivalence ratio, and temperature. 
However, there is a limiting equivalence ratio, for a given hydrocarbon-air system, above 
which solid carbon can form imder any condition. 

The formation of ammonia was considered only for the hydrogen-air system since 
the formation of solid carbon precluded the formation of ammonia for the combustion 
products of all the hydrocarbon-air systems studied. The limiting pressure for the for- 
mation of a given amount of ammonia depends upon the relative proportion of nitrogen N 2 
and hydrogen H 2 , equivalence ratio, and temperature. 

DERIVATION OF GAS MODEL 

Overall Stoichiometry and Gas Model Constants 

In order to apply the simplified gas model to the combustion of any hydrocarbon 
fuel, a generalized statement concerning the atomic composition of the fuel is required. 
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Consider the following stoichiometric reactions: 


CH + I O2 = CO2 + i H2O 
CH 2 + I O2 = CO2 + 1 H2O 
CH3 + 1 O2 = CO2 + I H2O 
CH4 + I O2 = CO2 + I H2O 

In general terms, 

CHn + O2 = CO2 + I H2O 

where in a fuel molecule 

No. atoms H 

n = 

No. atoms C 

Thus, in the general system, 

No. moles CHj^ 4rgq 
No. moles O 2 n + 4 

where r^q is the equivalence ratio. 

The initial composition of air incorporated in this gas model is the same as that of 
reference 3; that is, combustion is assumed to have occurred with dry air of the following 
composition by volume: 


O 2 20.9495 percent 

N 2 78.0881 percent 

Ar 0.9624 percent 

The following constants and parameters define, in general terms, the initial atomic pro- 
portions of the elements in the gas mixture: 

Q ^ No. atoms Ar ^ q. 00616227 
^ No. atoms N 

D 2 = No. atoms N ^ 3 727445 
No. atoms O 
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Do = No. atoms C _ 

^ No. atoms O n + 4 

£) _ No. atoms H _ 

^ No. atoms O n + 4 


and a convenient parameter 


Define r^.jj = i; then 


D = 1 + D2^1 + 2Di^ 


P3 - ^^eq^ch 


1 +4r 


ch 


2rpo 

°4 = 1 — 

^ 1 + 4r 


ch 


For the special system of hydrogen-air combustion, 


2H2 +02 = 2 H 2 O 

No. moles H 2 
No. moles O 2 

^ch " ® O 3 = 0 D 4 = 2rgq 


The stoichiometric fuel -air ratio is given by 


fg = 0.028931 


'1.008 + 12.01rgij) 

1 + 4rch 


It is of value to represent the equilibrium constants for the assumed reactions in a 
convenient manner. The treatment of chemical equilibria in reference 2 leads to the fol- 
lowing equilibrium expression for a general reaction: 


aA + bB = cC + dD 


The equilibrium expression is 




V-k 


X' 


B 
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where 


K = exp 





Figure 2 shows the variation of logjQ K with 1/T for the chemical reactions consid- 

Gt 

ered in this gas model. The values of Gibbs free energy ^ for the species were taken 
from the data tabulations in references 3 and 4. Figiure 2 indicates that the equilibrium 
constants can be represented as linear fimctions. 

With the initial proportions of the gas mixture elements defined and a convenient 
representation of the equilibrium constants determined, the governing relations for chem- 
ical composition can be formulated. Referring to figure 1 will aid the reader in the fol- 
lowing derivations. 


Hydrocarbon -Air Gas Model With Dissociation 
The following gas model applies to the region T > rgq > 0: 

Chemical species assumed 

H 2 O, CO 2 , O 2 , H 2 , N 2 , Ar, O, H, CO, OH 
Chemical reactions 
2H + O = H 2 O 
2H = H2 
20 = O 2 
O + CO = CO 2 
H + O = OH 

Mass conservation equations 
10 

(1) J Xj = 1 
J=1 

( 2 ) = T>i(2Xj^2) 

(3) 2Xn2 = D2(2Xqq2 + Xqo + ^HgO + 2 Xq^ + Xq + Xqh) 

Xc02 ^ ” ^3 (^Xc02 XCO + XH 2 O ^^02 Xq + Xqh) 

(5) Xjj + 2Xjj^ + 2 Xjj^q + Xqjj = D4^2Xq02 + X^q + + 2Xog + Xq + Xqjj^ 
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Equilibrium expressions 


(6) Xjj20 = p2KiXoX^ 


( 7 ) Xh2 = PK2x| 


(8) Xo2 = pK 3 X^ 


( 9 ) Xco2 = PVqXco 


(10) Xoh = PK 5 XoXh 


Solving for two equations and two unknowns in Xq and Xjj gives 

« = (a22Xo + A2iXo+ A2o)Xh + ^Ai2Xo 

+ Ai iXq + AiojXji + ^Ao 3 Xo + A02X0 + AqiXq + Aooj = 0 

(b22Xo + B21Xo + B2o)x^ + (bi2X^ ^ 

■ Bl iXq + Biojxn - ^Bq 3 X^ + B02X0 + BqiXo) = 0 

where 


A22 = p 3 KiK 4 (l - 2D3 + d) 

B22 = p 3 XiK 4 ^(l- 2 D 3 )-D 4 ] 

A21 = p 2 ^i(l + D) + 2K2K4(1 - 2D3) 

B2I = P 2 ^l[ 2 (l - D3) - D4] + 2K2K4 (i - 203)^ 

A20 = 2pK2(l - D3) 

B20 = 2pK2(^l - 03'^ 

Ai2 = P2X4K5(1 - 2D3 + d) 

B 12 = p2K4K5(l - 2D3 - D4) 

All = p^5(l + D) + 2X4(1 - 2D3) 

Bll = P^ 5 (l - B3 - D4) + K4 (i - 203) 

Aio = 2(1 - D3) 

Bio = 1 - D3 

Ao 3 = 2Dp2x3K4 

Bq 3 = 2 D 4 p 2 K 3 K 4 

Aq2 = p[2K3(d + D3) + K4(1 - 2D3 + d) 

Bq2 = B4p(2K3 + K4) 

Aoi = 1 + D - 2pK4(l - 2D3) 

BqI = D4 

Aoo = -2(1 - D3) 


Setting D3 and K4 equal to zero reduces the system to identically the hydrogen-air 

system. 
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The expressions for a and ^ are satisfied simultaneously by use of the Newton 
iteration scheme. Newton’s iteration scheme for determining Xq and Xjj can be 
written as 

+ 1) = Xo(L) - 

d^ 

where Xq(L) is the Lth approximation of Xq and the derivative of O! (Xo,Xh) with 
respect to Xq is 




da _ 

da 

da 


dXQ 

8Xq 

8Xh 



With Xq(L + 1) determined, Xjj(L + 1) can be calculated by applying the binomial 
theorem to the expression for i3. Although there may be a number of pairs of values 
^Xq,Xjj) which satisfy O! = 0, there is only one pair which satisfies the condition 
0 i Xj i 1. The initial value of Xq(L) is taken to be 0 and the iteration scheme is 
used luitil Xq is determined to five significant figures. 


For rgq > 1 the iteration function a is insensitive to changes in Xq. As a 
result, Xjj is used as the iteration variable in this region of computation. Unfortu- 
nately, if Xjj is assumed, the expression for /3 is cubic in Xq and an additional 
iteration is required. For fuel-rich hydrogen-air mixtures, Bq 3 = 0, ^ is quadratic in 
Xq, and the single-level iteration is maintained. However, for fuel-rich mixtures con- 
taining carbon, a two-level iteration is necessary for solution. 


Air Gas Model With Dissociation 

The following gas model applies to the region ^ ^ ^cut-off’ ^eq 
Chemical species assumed 
02, N2, At, O, NO, N 
Chemical reactions 
20 = O2 
2N = N2 
N + O = NO 
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Mass conservation equations 


6 

(1) ^ Xj = 1 

J^l 

(2) + ^NO ^n) 

(3) 2Xj^2 ■*■ ^NO “ ^2^2Xq2 + Xq + Xj^q) 

Equilibrium expressions 

(4) Xo2 = PK3X^ 

(5) Xnj^pKjX^ 

(6) Xj^Q = pK<jXqXj,j 

Solving for two equations and two unknowns in Xq and Xj^ yields 

a = A2 oX^ + (AuXq + Aio)Xn + (Ao2Xq + AoiXq + Aqo) = 0 

p = B20X§r + (BuXq + Biq)Xn - (bo2Xq + BoiXo) = 0 

where 

A 20 = (1 + 2Di)pKe B 20 = 2pK6 

All = (1 + Di)pK7 Bji = (1 - D2)pK7 

Aio = 1 + Di Bio = 1 


Aq 2 = PKs 


Bq 2 = 2 D 2 PK 3 


Aoi = 1 


BqI = E>2 


Aqo = 

The simultaneous solution of a and /3 is subject to the same iteration procedure used 
for the hydrocarbon-air gas model. 
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Fuel -Rich Gas Model With Solid Carbon Formation 
The following gas model applies to the region T ^ T(.ut-off> ^eq ^ 

Chemical species assumed 

H2O, CO2, H2, N2, At, CO 
Chemical reaction 

H2 + CO2 = CO + H2O 

Mass conservation equations 
6 

(1) J Xj = 1 

J=1 

(2) Xa = Di(2Xn2) 

(3) 2 Xj ^2 = D2^2 Xco 2 + Xco + XH20) 

^C02 + ^CO = 1^3(2X002 + ^CO + ^Hgo) 

(5) 2Xjj2 + 2 XJJ 2 O = E>4(^2Xco 2 + Xco + ^H2o) 

Equilibrium expression 

(6) Xg^o^CO = *^8^H2^C02 

Solving for Xqq yields 

AAXqq + BBXcq + CC = 0 

where 

AA = (1 - K8)(d + D4 + 2D3 - 1)^ 

BB = (d + D4 + 2D3 - 1) 2(1 - 2D3J + K8(D4 + 6D3 - 2) 

CC = 203X3(2 - D4 - 4D3) 

With the chemical composition in this region defined, the criteria for solid carbon 
formation can be formulated. Consider the reaction 

CO2 + C(Solid) = 2CO 

and equilibrium expression 

XcOo 
Plim = ^9 ^ 

■^CO 
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If p > solid carbon will form. There is, however, an equivalence ratio for the 

combustion products of a given hydrocarbon (r^h)* where solid carbon can form even in 
a near vacuum. Consider the following equilibrium expression: 

Pllm^CO “ *^*002 

(Plim-O) (XC02-°) 

Setting Xco 2 = 0 in the system for the fuel -rich gas model with solid carbon formation 
gives 

_l+4rch 
req,lim - 2rch 

Nondissociating Gas Model 
The following gas model applies to the region T ^ 

Chemical species assumed 

H 2 O, CO 2 , O 2 , H 2 , N 2 , Ar 
Mass conservation equations 
6 

(1) y Xj=l 
J^l 

(2) Xat = Di(2Xn2) 

(3) 2X^2 = D2^2 Xc 02 + X1J2O + 2 Xq 2^ 

(4) Xco2 = ^3(2X002 XH2O + 2X02) 

(5) 2 Xh 2 + 2XH2O = 04(2X002 + XH2O + 2X02) 

In each of the regions r^q = 0, 0 < rgq < 1, and rgq = 1, one or more of the species 

can be eliminated; thus, for rgq = 0, 

XH2 = 0 
Xco2 = 0 
XH 2 O = 0 
Xj = Constants 
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for 0 < r^q < 1, 


Xh2 = 0 

and for r^q = 1, 

Xh2 = 0 
Xo2 = 0 

Xj - Xj(rch) 

For the hydrogen-air system ^r^j^ = 0^ with r^q > 1, 

Xo2 = 0 

XC02 = ° 

“ ^j(^eq) 

With the chemical composition in this region, the condition for the formation of a given 
amount of ammonia can be determined. Consider the reaction 

N2 + 3 H 2 = 2NH3 
and equilibrium expression 

r2 


^ . ^NH3 

^10-^H2^N2 


(XnH 3 > 0.01, p > pj.^) 


Thermodynamic Properties of Gas Mixtures 

The previous derivations have outlined the solution for gas mixture composition. 
With the chemical composition defined and a tabulation of the thermodynamic properties 
of the pure constituents provided, the thermodynamic properties of a gas mixture can be 
calculated by the following equations (ref. 3): 
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where and are the thermodynamic properties of the pure constituents at refer- 

rC K 

ence pressure (1 atm). These thermodynamic expressions inherently assume that the 
pure constituents obey the perfect gas law. 

The standard reference state of the elements Ar, N, O, H, C taken from reference 3 
are as follows: 

Ar as Ar 

N as N 2 

O as O 2 

H as H 2 

C as C(Solid) 

By definition, the energy content of these elements in their standard reference states 
(T = 0° R) is zero. Equations for the computation of fuel enthalpies, consistent with this 
thermodynamic basis, are presented in reference 4. 

A listing of the computer program to calculate the composition and thermodynamic 
properties of this gas model is presented and discussed in the appendix. Computational 
time for obtaining the chemical composition and thermodynamic properties with this pro- 
gram has been estimated at 6000 cases/min on the IBM 7094 data processing system at 
the Langley Research Center. 


RESULTS AND DISCUSSION 
Comparison of Gas Models 

The computer program for the simplified combustion gas model has been used to 
calculate the thermodynamic properties of the stoichiometric kerosene -air and hydrogen- 
air systems. The purpose of this calculation is to compare the results of the simplified 
gas model with the more extensive treatments of references 3 and 4. Since the thermo- 
dynamic properties of the pure constituents used in the simplified model were taken from 
these references, the differences between the results of the simplified model and the ref- 
erence models are due to the formation of the neglected chemical species N, NH, NH 3 , 
and NO. 

Figure 3 is a Mollier diagram for the stoichiometric products of combustion of 
kerosene with air (r^jj = 0.5). The solid-line curves are values plotted from the tabulated 
data of reference 4. The dash-line curves are values obtained from the computer pro- 
gram for the simplified model. Unfortvmately, the data presented in reference 4 are lim- 
ited to temperatures below 5000° R (2778° K) . Since the formation of species including 
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atomic nitrogen ^N, NH, and NH 3 ) is more important at temperatures above 5000° PI, the 
good agreement between gas models is not surprising. 

Figure 4 is a Mollier diagram for the stoichiometric products of combustion of 
hydrogen with air ^r^j^ = oj. The solid-line curves are values plotted from the tabulated 
data of reference 3. The dash-line curves are values obtained from the computer pro- 
gram for the simplified model. The thermodynamic properties tabulated in reference 3 
cover temperatures as high as 10 000° R (5560° K). The substantial disagreement above 
7000° R (3889° K) at low pressures is due to the formation of species containing atomic 
nitrogen. An interesting result is that the agreement between lines of constant pressure 
is much better than the agreement between lines of constant temperature. This result is 
fortunate since, for thermodynamic processes such as isentropic nozzle expansion, tem- 
perature is not used directly to calculate performance. 

The simultaneous conditions of high temperatures and low pressures, the area of 
substantial disagreement between the simplified gas model and the reference gas models, 
are beyond the realm of operation of typical hypersonic engine cycles. For example, for 
subsonic combustion ramjets, combustion pressures and temperatures are high since the 
air stream is brought nearly to stagnation conditions before combustion. For supersonic 
combustion ramjets, combustion pressures and temperatures are low since a large por- 
tion of the total air enthalpy remains in the form of kinetic energy. Hence, the conditions 
of high temperatures and low pressures do not occvir simultaneously for engine cycles 
currently considered feasible for hypersonic flight. 

Ramjet Performance 

In applying the simplified gas model to hypersonic engine calculations, it is of inter- 
est to determine the effect of the use of the simplified gas model on the calculated per- 
formance of an idealized subsonic combustion ramjet. 

Since the purpose of this calculation is to show the effect of small differences in 
thermodynamic properties on ramjet performance, the following simplifying assumptions 
were made: 

(1) Free-stream Mach number chosen along a constant dynamic-pressure path 
(q = 1500 lb/ft2 (718 200 N/m2)) 

(2) Airstream decelerated to stagnation conditions with total pressure recovery 
degraded to 10 atm (due to internal duct pressure limitations at high Mach numbers) 

(3) Completely mixed stoichiometric hydrogen-air combustion 

(4) No pressure losses during combustion ^pj. nozzle " atm^ 


16 



(5) Enthalpy of injected hydrogen equals zero 

Molecular hydrogen is a reference element; thus, = 0 3^t T = 0° R. 

(6) Combustion gas isentropically expanded to free-stream static conditions by 
using so-called "shifting equilibrium" 

The 1962 Standard Atmosphere (ref. 7) yields free-stream static conditions and with 
assumption (1) gives the variation of altitude with Mach number shown in figure 5. For 
steady adiabatic flow, the total energy in the airstream is given by 

^t,air " ^«>,air 


For an adiabatic combustion process, 

„ _ ^t,air hydrogen 

^t, products “ 1 + f 

With two properties pj-Q^ju^ts Pt, nozzle) combustion gas defined, the 

Mollier diagram can be entered. Figure 6 is a schematic of a Mollier diagram showing 
lines of constant pressure for the simplified gas model and model of reference 3. Since 
the simplified gas model neglects the dissociated nitrogen species, the entropy level is 
slightly less than that obtained from reference 3. In expanding the combustion gas to 
free-stream static pressure, the gas enters the region of almost exact agreement between 
gas models. Consequently, the expansion using the simplified gas model gives a slightly 
larger value of AH = - Hg than that from reference 3. However, the effect of the 

difference in AH for the two calculations is reduced by the fact that 

Vj == ^Ht - Hs 

The performance parameters are given by the following expressions: 


T _ ^net 

®P w^gf 

The variation of specific impulse with Mach number for this idealized ramjet is presented 
in figure 7. The solid-line ctirve was calculated by using the thermodynamic properties 
of reference 3. The dash-line curve was calculated by using the computer program for 
the simplified gas model. The maximum deviation is 1 percent at a Mach number of 12. 
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Solid Carbon and Ammonia Formation 

The formation of solid carbon extracts useful energy from the combustion gas. If 
this phenomenon occurs during a nozzle expansion process, the useful energy absorbed in 
forming solid carbon is not recoverable and a loss in performance results. 

The derivation of the nondissociating ^T < fuel -rich gas model required 

one chemical reaction in order to include the formation of carbon monoxide CO. The 
equilibrium expression resulting from the required reaction is independent of pressure. 
Consequently, the composition of the gas in this region is independent of pressure. This 
fact simplifies the treatment of solid carbon formation, since the limiting pressure, above 
which solid carbon can form, is related to temperature and composition only. 

The computer program for the simplified gas model indicates the formation of solid 
carbon by an error statement. When the limiting pressure is exceeded, the error state- 
ment is printed, but the program computes the thermodynamic properties as though solid 
carbon had not formed. When the limiting equivalence ratio is exceeded, the error state- 
ment is printed and no calculations are made. The limiting pressure for solid carbon 
formation is readily extracted from the thermodynamic calculations. Figures 8 to 10 
present the variation of limiting pressure with equivalence ratio and temperature for 
methane -air = 0.25), kerosene -air (r^i^ = 0.5), and benzene -air (r^j^ = 1) combustion 
products, respectively. For all systems, as 

^eq ^eq,lim> Plim ® since — 0 

and as 

i*eq - h Plim •* ^CO ^ 0 

The formation of ammonia is considered only for the hydrogen-air, nondissociating, 
fuel-rich system. The composition of the gas, under these conditions, depends upon 
equivalence ratio only. Therefore, the limiting pressure, above which the mole fraction 
of ammonia is greater than some arbitrary amoimt, is related to temperature and 
composition. 

The computer program for the simplified gas model indicates the formation of 
ammonia > O.Ol) by an error statement. When the limiting pressure is exceeded, 

the error statement is printed, but the program computes the thermodynamic properties 
as though ammonia had not formed. Figiue 11 presents the variation of limiting pressure 
with equivalence ratio and temperature for hydrogen-air combustion products. For all 
temperatures, as 

^eq ^ Plim ^ ° 

and as 

^eq ^ 1. Plim ^ °° since ^ 0 
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CONCLUDING REMARKS 


A simplified equilibrium hydrocarbon-air combustion gas model, for use in engine 
cycle computer programs, has been presented. The generalized gas model includes the 
effects of dissociation and reduces to the special systems of hydrogen-air combustion 
products, as well as dissociating air. The associated computer program can be readily 
incorporated as a subroutine in a general engine cycle computer program or used alone 
to generate Mollier diagrams. 

With the exception of pure air, this gas model neglects the formation of chemical 
species containing atomic nitrogen; this assumption allows a considerable simplification 
of the solution for chemical composition. The effect of this assumption on the thermo- 
dynamic properties of stoichiometric kerosene-air and hydrogen-air combustion products 
is presented. The importance of this assumption is shown in terms of the performance 
of an idealized subsonic combustion ramjet. Good agreement between the simplified 
model and more comprehensive treatments is obtained in the range of temperatures 
applicable to hypersonic engine cycles. 

The computer program has been used to calculate the limiting pressure for solid 
carbon and ammonia formation in fuel -rich gas mixtures. The results of this calculation 
are presented as a fimction of temperature and equivalence ratio for the combustion prod- 
ucts of various hydrocarbon fuels. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 7, 1968, 

722-03-00-02-23. 
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APPENDIX 


COMPUTER PROGRAM FOR CALCULATION OF GAS COMPOSITION 
AND THERMODYNAMIC PROPERTIES OF 
SIMPLIFIED GAS MODEL 

The calculation procedure for determining gas composition and thermodynamic 
properties for the simplified gas model has been programed for and used with the 
IBM 7094 data processing system at the Langley Research Center. A printout of the 
program is presented in this appendix. The program is written in FORTRAN IV lan- 
guage (ref. 8). The symbols used for the program are as follows: 


RCH 

rch 

X(1) 

^H20 

REQ 

^eq 

X(2) 

Xc02 

TEMP 

T 

X(3) 

X02 

P 

P 

X(4) 

Xh2 

M 

M 

X(5) 

Xn2 

H 

H 

X(6) 

Xa 

S 

S 

X(7) 

Xo 

MT(J) 

Mj 

X(8) 

Xr 

TT(I) 

Tl 

X(9) 

Xco 

HT(I,J) 

H“j 

X(10) 

XoH 

ST(I,J) 


X(ll) 

Xno 



X(12) 

Xn 


Input 

The input is read into the IBM 7094 data processing system by the FORTRAN 
statement 

READ (5,100) RCH, REQ, TEMP, P 
100 FORMAT (4E 16.8) 

For example, an input car would be 

Col.- 1 17 23 49 

+0.00000000E+00 +0.10000000E+01 +0.50000000E+04 +0.10000000E+00 
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Output 

The outputs for this program are the mole fractions of each constituent and the 
thermodynamic properties of the gas mixture. For example, output for the input example 
would be as follows: 

RCH=O.OOOOOOOOE-38 REQ=0.10000000E+01 TEMP=0. 50000000 E 04 R P=0.100000000E00ATM 


XH20 

0.18315569E 

00 

XC02 

O.OOOOOOOOE-38 

X02 

0.27084661E-01 

XH2 

0.71114589E-01 

XN2 

0.57613321E 

00 

XA 

0.71005768E-02 

XO 

0.25695777E-01 

XH 

0.63610586E-01 

XCO 

O.OOOOOOOOE- 

-38 

XOH 

0.46109576E-01 

XNO 

O.OOOOOOOOE-38 

XN 

O.OOOOOOOOE -38 


M= 0.2 199386 5E 02 H=0.14316146E 04BTU/LBM S=0.30918686E OIBTU/LBM-R 


Built-In Data 


The thermodynamic properties of the pure constituents -^1 are taken from ref- 


erences 3 and 4 and built in the program as two-dimensional arrays named HT(I,J) and 
ST(I, J), where I is the index on temperature and J is the index on constituents. The 
thermodynamic properties are built in for the discrete temperatures given in TT(I). 

The enthalpy array HT(I, J) is interpolated linearly, whereas the entropy array ST(I,J) 
is interpolated logarithmically. 


Complete Program 

The complete program including comments is reproduced in the following pages. 
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= H-fX(J)*(hT(I-l,J) + (HT(I,J)-HT(I-l«J))'*-DELH) 
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Figure I.- Description of gas model. 
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Figure 2.- Variation of equilibrium constants with temperature for gas model reactions. 
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Enfholpy, H, Btu/lbm 



Entropy, S,J/kg-®K 


Figure 3.- Mollier diagram for kerosene-air combustion products (r^h ~ 0-5). tgq = 1.0. 
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Figure 4.- Mollier diagram for hydrogen-air combustion products = o). 
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Figure 5.- Assumed flight profile for ramjet performance calculation. 




Ref. 3 

Simplified model 



oas ‘ dsj 


42 


Figure 7.- Effect of gas model on typical ramjet performance. Hydrogen-air combustion; r^q - 1.0; q = 1500 ibf/ft^ (718 200 N/m^); 
p. - 10 atm; isentropic fully expanded nozzle; 1962 U.S. Standard Atmosphere. 
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Figure 8.- Limiting pressure for solid carbon formation for methane-air combustion products (r^^ = 0.25). rgg |jj^ = 4.0. 
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Figure 9.- Limiting pressure for solid carbon formation for kerosene-air combustion products (r^h = 0.5). r^g 
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Equivalence ratio, r^q 


Figure 10.- Limiting pressure for solid carbon formation for benzene-air combustion products (r^h = l-O). r^q |jm = 2.5. 
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Figure 11.- Limiting pressure for ammonia formation = 0.01^ for hydrogen-air combustion products (r^-^ = 0). 
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